Matrix rearrangement approach for the entangling power with mixed qudit systems 
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We extend the former matrix rearrangement approach of the entangling power to the general 
cases, without the requirement of the same dimensions of the subsystems. The entangling power 
of a unitary operator is completely determined by its realignment and partial transposition. As 
applications, we calculate the entangling power for the Ising interaction and the isotropic Heisenberg 
interaction in the mixed qudit system. 



I. INTRODUCTION 

Entanglement is one of the most important quantum correlations among subsystems. It plays a key role in quantum 
information processing and is considered as a kind of resource for quantum computation and communication In 
recent decades, many efforts are being put into the quantification, generation and enhancement of the entanglement. 

From the aspect of generation of entanglement, what we concern is that given an operation, what is its capability of 
producing entanglement. Recently, the entanglement capability of quantum evolution and Hamiltonians is studied in 
Refs. @, H, 0, H, 0, 0, H, @ • Generally speaking, for unitary operators, we can let them act on pure product states, then 
observe the amount of the entanglement produced. It always depends on the input states. There have been two ways 
to obtain a input-state- independent quantity for measuring the unitary operators' capacity of entanglement creation. 
One is averaging over the product input state 0, [Tot [ill . [l2| , the other is taking the maximum of the entanglement 
produced over the product input states [H, • 

The entangling power belongs to the first kind. It is defined as [l(| 



e v {U) = E{U\j>i) <8> (1) 

where the average is over all the product states distributed according to some probability distribution. The linear 
entropy E is adopted as entanglement measurement for pure states of bipartite systems hereafter, which is defined as 

£(|*» = 1-TrM), (2) 

where p\ = Tr2(| 1 S , )(W|) is the reduced density matrix of the first subsystem. 

When we study the closed quantum systems, under Schrodinger equation, the time evolution is unitary. Especially, 
when the Hamiltonian is time- independent, the unitary time evolution operator is in the one-parameter group gener- 
ated by the Hamiltonian, i.e., U(t) — exp(— iHt). So we can also consider conveniently the time average of entangling 
power which is define as 

T 

e^(U(t)) = r lim - J e p (U(t))dt. (3) 
o 

This entangling power not only average over all initial product states, but also over all time range. 

Most recently, a matrix rearrangement approach to the entangling power was given in [l5| for subsystems with 
equal dimensions. It expresses the entangling power in terms of the realigned matrix [l6| and the partially transposed 
matrix [13]. These two matrix rearrangements are originally used to study the separability problem of the quantum 
mixed states, and lead to two useful criterion for entanglement fl6l. [T7L [l8l [l9| . 

On the other hand, because the entanglement plays an intriguing role in the kinematic process of the reducing of 
density operators [2(J, [HJ and the reduced dynamics (23, the entangling power has also been used to study some 
dynamical systems, especially in the chaotic systems [23l |24|. 

In this paper, first, we prove that the entangling power of an unitary operator U is completely determined by its 
realignment and partial transposition, without requirement of same dimensions of the subsystems. Then we apply 
it to calculate the entangling power of the Ising interaction and the isotropic Heisenberg interaction in hybrid qudit 
d\ x d,2 systems. The analytic results are obtained. We can see the exact relation between the time-average entangling 
power and dimensions of the Hilbert spaces of the subsystems. At last, a conclusion is given. 
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II. MATRIX REARRANGEMENT APPROACH FOR THE ENTANGLING POWER 

After averaging with respect to the uniform distribution of the non-entangled pure initial states, for a unitary operator 
on d\ x g? 2 systems, the entangling power is given by fioj 

e p (U) = l-C dl C d2 *a(E0, 

a=0,l (4) 

I a (U) = Tr(T 1+Q , 3+ct ) + {U® 2 , (T 1+a , 3+a )U® 2 T 13 ), 

where C^ 1 — d(d + 1), and {A, B) = Tr(A^B) is the Hilbert-Schmidt scalar product. Let Tjj = 1, . . . , 4) denote 
the transposition between the ith and the jth factor of H® 2 = {C dl ® C d2 ) (g> (C dl (g> C d2 ). For example, 

2i 3 (|*> ® |j) ® ® 10) = I*) ® b'> ® l») ® 10- ( 5 ) 
It is easy to check that Tr(Ti3) = did 2 , and Tr(T24) = d\da. Then we get another form of the entangling power 

e P (U) = T i 1: d 2 ~ g dl ,d2 E ( U ®*> (Ti +a , 3+a )U® 2 T 13 ), (6) 



a=0,l 



where 



d x d 2 + 1 
TdlA = (di + l)(da + l)' 

^ = did 2 (di + l)(da + 1) " (?) 
For the dxd systems (di = (fe), the entangling power can be expressed as [25[ 

gp = { dTT )2[E{u) + ^ (c/5l2) " E{Sl2)] > (8) 

where -E(O) is the entanglement of the operator O and can be expressed using the matrix realignment as E(0) = 
1 — Tr[(U R (U R )^) 2 ]/d 4 , so the entangling power has the following form [l5[ 

e P (U) =^±^ - ^^{Tr[(^(^) t ) 2 ] +Tr[(^(f/ Tl ) t ) 2 ]}, 0) 

where R is the matrix realignment [l6| and Ti the partial transposition with respect to the first subsystem (l7| . They 
are defined by 

(U R )ij,ki = {U)ik,jh 
(U Tl )ij,ki — {U)kj,u, 

respectively. 

We now show that the entangling power is determined by the realigned and partially transposed unitary operator 
without requiring the same dimensions d\ = d%. From Eq. we see that the terms related to the evolution operator 
are only (U® 2 , Ti +a ^ +a U® 2 Ti 3 ) . So the task left is to get the other expressions of them. A unitary operator can be 
written as 

U = Uiajpeij ® e a/ 3 

= A-ij,a/3 e ij ® S a p 

= Tja t ipeij®e a p, (11) 

where e%j = Repeated indices are summed in all cases. Here we use English letter subfixes for the first 

subsystem, and Greek letter ones for the second. Obviously, in the expression above, A = [/ fl ,and T — U Tl . Then we 
get 

U® 2 = k iha pk k i n $eij ® e a p ® e k i ® e 7< 5 
= Fja^pTij^seij ® e Q( 3 ® e k i ® e 7< 5, 
Ti 3 U® 2 Ti3 = A.ij, a pAki,jseki ® e Q/3 ® ® e 7<5 , 

T 2 aU® 2 T 1z = Tja^pTi^^seu ® e 7/3 ® ejy e Q<5 . (12) 
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Substituting Eq. (|12p into the related terms, and notice the normalization relation {eij,e k i) = SikSji, we get 

(U® 2 ,T 13 U® 2 T 13 ) = A^A^A^A^ = Tr[(AA^) 2 ], (13) 

and 

(U® 2 ,T 24 U® 2 T 13 ) = T* ja , ip T^ tkS T hiifi T jchkS = Tr[(ITt) 2 ]. (14) 

Substituting them into Eq. ([6]), note that A = U , and T = £/ Tl , we get the general expression of the entangling 
power 

e P (U) = T dlM - g dl , d2 {TrW R U Ri ) 2 ] + Tr[(U T ^) 2 }}. (15) 

That means for the systems with certain dimensions, the entangling power is proportional to the trace of the product 
of the matrices after realignment and partial transposition to the time evolution operator. In other words, we have 
shown that the entangling power is completely determined by the realigned and partially transposed operators. These 
two matrix manipulations are sufficient to study the entangling capability of an unitary operator. The only task to 
get the entangling power is to calculate these trace terms. 

For later uses, we provide a simple form of these trace terms. We can always decompose the evolution operator 
into the form: 

U = Y^<kA i %B i . (16) 

i 

After some algebra, we get 

Tr[(U R U Rt ) 2 } =^c i c*c fcC rTr[A iJ 4/]Tr[^ fe A ; t ]Tr[B l B ; t ]Tr[B fc B/], 

Tr[(£/ Tl C/ TlT ) 2 ] =J2ciC*c k c^Tr[A i A^A k A^]Tr{B i B l ^B k B^]. 

ijkl 

We can intuitively see the connection between the entangling power and the decomposition of the evolution operator 

here. If the evolution operator can be written as U = J7i<& U 2 , we have Tr[(U R U R *) 2 ] = d\d\, Tr[(£/ Tl t/ Tlt ) 2 ] = d\da, 
then from Eq. (|15p . the entangling power vanishes. Equation (|17p will be used when calculating the entangling power 
of the Heisenberg interaction in Sec. IV. 



III. APPLICATION Is ISING INTERACTION 



We study the entangling power of the time evolution of the Ising model, whose Hamiltonian reads H = gSl ® S%- 
The unitary time evolution operator generated by it is 



si s 2 



U= E e Wm ^\m im2 )(m im2 \, (18) 



m 1 = — s 1 m 2 = -s 2 



where 8 = —gt, Si is the length of spin i and |mj) is the eigenstate of Sf. ft = 1 is assumed for all cases. 
Making use of the definition of realignment and partial transposition, we get 

U Tl = U, (19) 
U R = e i0mim2 \m 1 m 1 )(m 2 m 2 \. (20) 

7711 771-2 

Then 

U R (U R )t= (Xy e(mi ~™ l)m2 )K TO i}(™i n il- ( 21 ) 

mim.2 Tii 
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Note that the range of the indices in the above expression is all symmetric, so we can substitute index m with — m. 
Then after trace, we get 



Tt[(U r (U r )^) 2 } 

si s 2 

= E ( E 

m\,n\— — s\ m2— — S2 
di-1 



a i8(m 1 +n 1 )m 2 \2 



(22) 



h<% + 2(di-M) — 



M=l 



1 - cos(d 2 M9) 
cos(A/<9) ' 



and 



Tr[(C/ Tl (i7 Tl ) t ) 2 ] = Tr[(UU f ) 2 } = did 2 , (23) 
where di = 2s t + 1 is the dimensions of the Hilbert space of spin i. Substituting Eqs. (|22| and ([23]) to (fl~5"j) leads to 

, m (ji - 1)^2 1 1 - cos(d 2 M6>) 

6p(C/) = (di + !)(* + !) " d ld2 (d 1 + l)(d 2 + l) E 2(* ~ M ) ! ~— • ( 24 ) 



(di + l)(d2 + l) ^2(^ + 1)^2 + 1)^ 

For the cases of the qubit-qudit system (di = 2), we get 

d 2 1 - cos^f?) 



cos(AJ0) 



e P (U) 



3(d 2 + l) 3d 2 (d 2 + 1)(1 -cosfl)' 
And for the cases of two-qubit system (di = d 2 = 2), we get the same result as [HI 

e p (lO = ^sin 2 (0/2). 



(25) 



(26) 




FIG. 1: The entangling power of the Ising model. In the left figure we let di = 2. 

From Fig. (fTJ, we may see that when the dimensions of the subsystems increase, the entangling power increases on 
average, and the time oscillations become more complex when dimensions increase. At 9 = 7r, the entangling power 
reaches its extremum, but often it is a minimum other than maximum. The maxima of entangling power mean large 
entanglement generation, and these points correspond to some useful quantum gates. So, the study of entangling 
power helps us to construct quantum gates in quantum computation with hybrid qudits. 

We further consider the time average of entangling power defined by Eq. ([3]) . Because 



Z7T 

— / Ti[{U R {U R )^) 2 ]d9 
2n J 

o 

~y ^ y ^ I e i9(rn 1 +ni)(m2+n 2 ) 



m!,ni = — sx m 2 ,n 2 =— s 2 . 

d 1 d 2 (d 1 +d 2 -l), (27) 
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then after substituting it to the definition of the time-average entangling power Eq. ([3]) , we get 

^) = (1-^-t)(1-^-t)- (28) 
d\ + 1 d 2 + 1 

It increases as the dimensions of the subsystems increases, which can also be seen from Fig. (p} intuitively. 

IV. APPLICATION II: HEISENBERG INTERACTION 

We now consider a general SU(2)-invariant Hamiltonian. According to the Schur's lemma, it can be written as 

S2+S1 



H= E ^ ( 29 ) 



where P n is projection operator of the total spin-n subspaces and E n is the eigenenergy. Without loss of generality, 
s i 5= s 2 is assumed here. The time evolution operator generated by it is given by 



S2+S1 



u= a « p «' ( 3 °) 



where a n — e ltEn . 

On the other hand, the spin projection operators can be expressed in term of the SU(2) operators (Si • S2)™ as 

dssa 

Pn = n \ ; (3D 

K—S2—S1 

where 

A fe = i [k(k + 1) - si(ai + 1) - s 2 (s 2 + 1)] • (32) 
For instance, in the case of si = 1/2, we get 

d - 1 - 4Si • S 2 



P, 



S2-1/2 



2d 



d+l + 4S 1 -S 3 

i^+i/a - Yd ' (33) 



where d — 2s 2 + 1 is the dimensions of the Hilbert space of the second spin. 
Substituting Eq. (f3Tj) into (|30p. we formally get the following useful expression 



t/ = ^/3„(S 1 -S 2 )™ ! (34) 



n=0 



where {3 n is the corresponding parameter. Substituting Si • S2 = J2i=x v Z S\ ® S 2 into the above expression and 
expand it, then we get the decomposition form as in Eq. (|16p . 

For an isotropic Heisenberg interaction bipartite model, the Hamiltonian reads H = Si • S 2 . In the case of Si = 1/2, 
due to Eq. ((33]), we can write the Hamiltonian as 

H = E P S2 _ 1/2 + £iP S2+1/2 , (35) 

where Eq = —(d+ l)/4, E\ = (d— l)/4 are the eigenvalues of the Hamiltonian. And we can write the time evolution 
operator as follows 

U = ftl + ftSi-Sa, 

/3b = [(d-l)e- itEo + (d+l)e~ ztEl ]/{2d), (36) 
/3i = 2{e- UEl -e- UEo )/d, 
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where I is the identity operator. Making use of Eq. (fT7|) , and noticing that 

Tr(S^) = Sa £ fc 2 = ~ X K ( 37 ) 

k— — s 

Ti(S l ) = 0, (38) 

where i € {x, y, z}, we get 

Tr[(U R U^) 2 ] = 4d 2 |/3 | 4 + -^{d 2 - l) 2 d 2 \(3i\\ (39) 

On the other hand, for the partial transpose operation, we can use the partial time reversal [28[ instead, which is 
defined by 

M T1 = e-™ s *M Tl e™ s K (40) 

We have 

Tr(([/ Tl C/ Tlt ) 2 ) = Tr(([/ T1 C/ Tlt ) 2 ). (41) 

Then we get 

C/ ri = /3 /-/3 1 Si-S 2 =7oP +7i^i, (42) 
7o = fa-Eofa, j 1 =p -E 1 (3 1 , 

where Pq = P S2 -i/2 and -Pi = P S2 +i/2 are the projection operators of the subspaces with the eigenenergy Eq and E\, 
respectively. Hence we obtain 

Tr[(U T ^ j ) 2 ] = (d 1)| 70 | 4 + (d+ l)|7i| 4 - (43) 
At last, we get the entangling power 

6p(C/) = A -^dT- [M " ~ /W^ 2 (|)]sin 2 (|), (44) 

where / (d) = 2 (6 — d + d 3 ) . In the special case of d = 2 , which means a two-qubit system, the above equation reduces 
to 

ep (C/) = isin 2 (i), (45) 
6 

which is consistent with the entangling power of a SWAP operator = 2Si ■ S2 + 1/2 [15| . And the time-average 
entangling power is 

(d- l){d 3 + d-6) 
3d 4 

It increases when the dimensions of the subsystems increases, as the case of the Ising interaction. 

We also make numerical calculations which are shown in Fig. @. We can see that the period of the entangling 
power is smaller when d is greater. This is due to the propertis of the time evolution operator. Because of the SU(2) 
invariance, the Hamiltonian describing the isotropic Heisenberg interaction between a qubit and a qudit system has 
only two different eigenvalues, and the energy gap becomes greater as the dimensions of the Hilbcrt space of the qudit 
subsystem increases, then the period of the time evolution becomes smaller. 



~ P (U) = V " ^7 (46) 



V. CONCLUSION 



In conclusion, we have extended the matrix rearrangement approach for the entangling power to the general cases of 
the hybrid qudit (di x d2) systems. This approach supplies a convenient way to get the entangling power, and we show 
that the exact solutions of the entangling power are obtained for some simple Hamiltonians. We also consider the 
effects of the dimensions of the subsystems, and find that the time average entangling power are monotone increasing 
functions with respect to the dimensions of the subsystems. 

Comparing the formula given by Zanardi et al. ,10|, the present one is simpler to apply. What we need to do is only 
to calculate the realigned matrix and the partially transposed matrix, and after making some simple traces, we can 
obtain the entangling power. The entangling power has been successfully used in the study of quantum chaos [291 ]. and 
we believe that the time-average entangling power introduced here is also useful for characterizing nonlinear behaviors 
of quantum systems. 



t 



FIG. 2: The entangling power of the isotropic Heisenberg model. 
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